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Abstract 



In this paper, convergent, multi-infinite, series solutions are derived for the homoclinic 
orbits of a canonical fourth-order ODE system, in both reversible and non-reversible cases. 
This ODE includes traveling- wave reductions of many important nonlinear PDEs or PDE 
systems, for which these analytical solutions would correspond to regular or localized 
pulses of the PDE. As such, the homoclinic solutions derived here are clearly topical, 
and they are shown to match closely to earlier results obtained by homoclinic numerical 
shooting. In addition, the results for the non-reversible case go beyond those that have 
been typically considered in analyses conducted within bifurcation-theoretic settings. 

We also comment on generalizing the treatment here to parameter regimes where 
solutions homoclinic to exponentially small periodic orbits are known to exist, as well as 
£C) ■ another possible extension placing the solutions derived here within the framework of a 

comprehensive categorization of ALL possible traveling-wave solutions, both smooth and 
non-smooth, for our governing ODE. 



1 Introduction 



In this paper, we initiate a novel analytic approach to the calculation of homoclinic orbits in 
both reversible and non-reversible ODE systems. 
Homoclinic orbits of dynamical systems have been widely treated in recent years by a variety 
of approaches. For instance, an early review integrating bifurcation theoretic and numerical 
approaches was given in pQ. Homoclinic orbits are important in applications for a variety of 
reasons. In the context of ODE systems, they are often anchors for the local dynamics in 
their vicinity. Under certain conditions, their existence may indicate the existence of chaos 
in their neighborhood (see [3] - [I] for instance). 
In a totally different setting, if the governing dynamical system is the traveling-wave ODE 
for a partial differential equation or equations, its homoclinic orbits correspond to the solitary 
wave or pulse solutions of the PDE(s), which have many important uses and applications in 
nonlinear wave propagation theory, nonlinear optics, and in various other settings (see [5] - 
[6] for instance). 

Given their importance, we will develop an analytic solution which will cover, or be pertinent 
to, diverse systems of interest in the relevant regime of parameters. Towards this end, we start 
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with an important, canonical ODE system which may, among other applications, result from 
traveling-wave reductions of a variety of physically-important nonlinear PDEs (henceforth 
written NLPDEs). As such, this ODE has served as an important workhorse in treatments 
of homoclinic orbits in reversible systems pQ. We shall also generalize our treatment and 
include non-reversible terms (nonlinearities) which have been considered only rarely, if at all, 
in earlier work. 

This generalized canonical ODE is introduced in Section 2, and followed in Section 3 by 
a review of some known results for its homoclinic solutions in four complementary regions 
spanning all of the parameter space. 

Of these regions, two contain countably infinite numbers, or a so-called plethora, of homo- 
clinic orbits. In one of these two regions, the orbits are homoclinic to the zero solution. We 
derive explicit analytic solutions for the homoclinic orbits in that region in Sections 4 and 
5 for the reversible and non-reversible cases respectively. If the ODE is a traveling-wave 
reduction of some PDE, these solutions would correspond to genuine or localized pulses of 
the PDE. 

In the second of the above regions containing an infinity of such orbits, the homoclinic 
solutions are known to be homoclinic or asymptotic to exponentially small periodic orbits. 
These correspond to delocalized solitary traveling wave solutions (with exponentially small 
tails) of any corresponding PDE(s). We do not treat the solutions in this second region in 
this paper, but defer them to a future report. 

The remainder of the paper considers our analytic homoclinic solutions, which have the form 
of multi-infinite series, in greater detail. So as not to affect the readability of the paper, we 
defer a convergence proof for our series solutions to Section 6. Before that, the solutions are 
numerically summed and plotted in Sections 4 and 5 - fortunately (and this is subsequently 
proven using an inductive argument as part of the convergence proof of Section 6), the later 
coefficients drop off sharply in magnitude (which is by no means guaranteed a priori for such 
convergent series). The plots in Section 4 are in very close agreement with earlier numerical 
results for the reversible case obtained by homoclinic shooting. As for the corresponding 
homoclinic orbits derived for the non-reversible case in Section 5, to our knowledge, the 
results there go beyond most which exist in the literature. 

For completeness, we should mention at this point that the only OTHER existing results 
(sec |5j-[6j for instance) for homoclinic orbits of our ODE, besides the numerical solutions 
mentioned above, are from the use of the bifurcation-tracking software AUTO. Although it 
would be straightforward to compare our series solutions to those as well, we do not consider 
them as they have been numerically benchmarked earlier against those obtained by homoclinic 
shooting, which we do compare our results against. 

Section 6 rounds out the paper with a proof of the convergence of our analytical multi- 
infinite series solutions for the homoclinic orbits, while Section 7 summarizes the results, 
and points out directions for further investigations. In particular, we make some comments 
on (non-trivially) extending our treatment to the second region mentioned above where the 
solutions are homoclinic to exponentially small periodic orbits. We also remark briefly on a 
forthcoming paper where the solutions and results here will be juxtaposed with, and placed in 
the framework of, a comprehensive categorization of ALL possible traveling-wave solutions, 
both smooth and non-smooth, for our governing ODE. 
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Figure 1: Linearization at the origin of (|2.ip . 

2 A canonical fourth-order equation and related physical PDEs 

Let us consider the equation: 

//// // . / // /// , . . 

u — bu + au = j (u,u ,u ,u ), (2.1) 

where a and b are real parameters, / is a nonlinear function, and the ' indicates the derivatives 
with respect to z. 

As one possible motivation, recall that substituting traveling wave solutions of the following 
form: 

u(x,t) = u(x — at) = u{z) (2-2) 
into the fifth-order KdV5 or FKdV equation: 

ut + 6uu x + u xxx + u xxxxx + 10uu xxx + 20u x u xx + 30u 2 u x = (2-3) 
as done, for instance, in [2], one obtains the following equation: 

//// // Q // /r) tj 

u + u - au + 3u 2 + Wuu + 5u 2 + 10n d = (2.4) 

which belongs to the class (|2.ip . Numerous other NLPDEs of importance in various applica- 
tion areas and having traveling wave-reduced equations of the form (I2.ip are discussed in [Tj 
and the references there. 



3 Possible solitary wave domains via bifurcation analysis 

In order to examine the structure of homoclinic solutions to (12. ip the first step is to review 
its linearization pQ, which is depicted in Fig{TJ 
In FigQ] there are four distinct regions, bounded by the codimension one curves Cj, i = 
0, . . . , 3, and corresponding to qualitatively different linear dynamics in each part. 
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Equation (|2,ip is invariant under the transformation z — > —z and is thus a reversible system. 
In this section, we shall use the theory of reversible systems to characterize the homoclinic 
orbits to the fixed point of (|2.ip in various regions of the (b, a) plane. These would correspond 
to pulses or solitary wave solutions of any NLPDE (such as those discussed above) for which 
(|2.ip is the traveling- wave ODE, . 

The linearized system corresponding to (12. ip 

mi it . , 

u — bu + au = (3-1) 

has a fixed point 

u = u = u = u = 0. (3-2) 

Solutions 4> = ke Xz satisfy the characteristic equation A 4 — 6A 2 + a = from which one may 
deduce that the structure of the eigenvalues is distinct in four different regions of the (6, a) 
plane. These are the regions shown in Figure [Hand labeled regions 1-4. 
In delineating the structure of the homoclinic orbit to fixed point (|3,2p in various parts of 
the (b, a) space, we shall first consider the bounding curves Co — C3 and their neighborhoods. 
Following this, we shall discuss the possible occurrence and multiplicities of homoclinic orbits 
to (13. 2p in each of regions 1 through 4: 

a. Near Co: This curve, on which the eigenvalues have the structure Ai_4 = 0, 0, ±A, and 
its vicinity have been considered in the context of reversible systems in [Z]-[H]- In this 
region, a standard analysis yields the normal form on the center manifold 

x\ = x 2 

x 2 = sign(/i)xi - -x\ 

where [i is an unfolding parameter [7J. For > 0, this yields a unique, symmetric 
homoclinic solution 

x x {t) = sech 2 (t/2) 

in the vicinity of Co- One may also show persistence of this homoclinic solution in the 
original system (2.1) for \x > [7]. 

b. Near C\: Near Ci, which corresponds to the eigenvalue structure Ai_4 = 0, 0, ±iuj, we 
will show by analysis of a four-dimensional normal form [7J in Section 4 that on the 
side of C\ corresponding to region 3 in Figure 1 there is a sech 2 homoclinic orbit. 

However, in region 3, where the eigenvalue structure is that of a saddle-center Ai_4 = 
±A,zw, the fixed point (j3.2j) is non-hyperbolic. It may be proved [9]-[TC)j that there 
are classes of homoclinic orbits in region 3 which are homoclinic to periodic orbits as 
z —> ±00. These are so-called delocalized solitary waves |11| . Depending on the form 
of the nonlinear term (the RHS of (2.1) in our case), these periodic orbits often have 
exponentially small amplitudes. Also, on isolated curves in region 3, the amplitude of 
these periodic solutions as z — > ±00 goes to zero, thus yielding truly localized solitary 
waves. These are known as embedded solitons [12J and will be investigated for our 
system (2.1) in future work. 
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c. Near C^- In this region, where Ai_4 = ±iu), derivation and analysis of a com- 
plicated normal form p3]-[T3] in Section 5 shows the possible occurrence of so-called 
envelope homoclinic solutions of the form sech(/ct)e^ 7e ^ (and with oscillating tails) in 
the so-called subcritical form. However, occurrence or persistence of these solutions in 
the full nonlinear system (2.1) is a non-trivial issue (and each system must be analyzed 
separately |14j-|16j). The persistence in (2.1) is considered in various ways in |14j-|16|. 
Indeed, there is no problem for solutions having one hump or peak. The only open 
problem is for nonsymmetric solutions. 

d. Near C3: There is no small-amplitude bifurcation on C3, on which Ai_4 = ±A,±A 
and the fixed point ()3.2[) remains hyperbolic. However, as we discuss below, there is a 
bifurcation across it causing the creation of an infinite multiplicity of homoclinic orbits. 

We turn next to each of the regions 1 to 4 in Figure 1 to discuss the possible occurrence and 
multiplicity of homoclinic orbits in each one. 

a. Region 1. In this region Ai_4 = ±\±iu and fixed point (13. 2j) is a saddle-focus. Using 
a Shil'nikov type analysis, one may show [17j.[18| for general reversible systems such 
as (2.1) that the existence of one symmetric homoclinic orbit implies the existence of 
an infinity of others. Hence, we expect our system (2.1) to admit an infinity of such 
symmetric A^-pulses for each N > 1. Here, a symmetric iV-pulse oscillates N times in 
phase-space for z £ (—00,00) (or, more technically, crosses a transversal section to the 
primary symmetric 1 pulse iV times). In the context of the Ostrovsky equation, these 
would be TV-peaked solitary waves, and we expect an infinite family for all N > 1 for 
parameters (p, q) in Region 1 of Figure 1. 

b. Region 2. In this region Ai_4 = ±Ai, ±A2 and fixed point (|3.2p is a hyperbolic saddle 
point. Thus, there is no a priori reason for multiplicity of homoclinic orbits in this 
region. However, depending on the actual form of the nonlinear term, a symmetric ho- 
moclinic orbit to (|3.2p may exist (see [I3],[I1]). Also, depending on further conditions 
|15j.[16j. a further orbit-flip bifurcation may cause complex dynamics in its neighbor- 
hood. In the context of our system (2.1), these issues will need further investigation to 
establish possible existence of solitary wave solutions in this region of (b, a) space. 

c. Region 3. The generic situation in this region has already been considered in the 
discussion above pertaining to the region near curve C\. As mentioned there, the 
structure and multiplicity of the delocalized solitons in Region 3, as well as the existence 
of embedded solitons on isolated curves, will be investigated for the Ostrovsky equation 
in future work. 

d. Region 4. In this region, Ai_4 = ±iwi, ±iw2 and (|3.2p is a focus. No homoclinic orbits 
are known to exist in general here, although complex dynamics may occur |17j.|18j. 
Special results exist for uj\ « u>2- 

Having outlined the possible families of orbits homoclinic to fixed point (I3.2p of (|3.ip . we 
proceed next to derive convergent series solutions for solutions homoclinic to this fixed point 
in Region 1 discussed above. 
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4 Reversible case: symmetric homoclinic solutions 



Here we shall derive solutions in the case where the nonlinear terms are chosen so that (|2.ip 
is reversible as discussed in the preceding section. However, the method of undetermined 
coefficients we employ for the analytic solutions proves robust, and hence, in the following 
section, the series solutions developed here are generalized to include nonlinear terms such 
that (]2.ip may be non-reversible, and hence to cases more general than may be treated within 
the qualitative, reversible systems bifurcation framework which has been discussed up to this 
point. 

In this section, we choose the nonlinear function / in (|2.ip in such a way that the equation 
is reversible. So as to include the traveling wave reductions of various, physically-relevant 
nonlinear PDEs discussed in pQ, we choose the nonlinearity to be of the form: 

. I II III. tl /O 2 T '"i / \ 

j(u,u,u,u ) = cuu + du + gu + hu . (4-1) 

The resulting equation: 

u — bu + au = cuu + du + gu + hu , (4-2) 

is reversible under the standard reversibility of classical mechanical systems: 

i u ut t a in 

z — » —z (u,u ,u ,u )^-{u,—u,u ,—u ). (4-3) 

Mathematically, this translates to solutions having even parity in z. 
The equation (|4.2j) has the following three fixed points: 



u = 0, u= y y ^ , MO (4.4) 

Notice that if h = (i.e. if there is no u 3 term) there are only two fixed points u = and 
u = a/g. 

Let us now proceed to construct in region 1 of Figure 1 Shilnikov-type homoclinic orbits to 
the origin. Suppose that for z > 0: 

oo 

u+(z) = ^2a k e kaz , (4.5) 

k=l 

where a < is an undetermined constant and a/-, k > 1 are, at the outset, arbitrary coeffi- 
cients. Substituting the series f|4.5j) into the reversible equation f|4.2j) yields: 

oo 

{(ka) 4 - b(ka) 2 + a) a k e kaz = F x + F 2 , (4.6) 

k=l 

with: 



Fl = J2Fi k \ F 2 = f2H k) (4-7) 



f), /•:, Y, 

k=2 fc=3 

where 
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{0 for k = 1 

(4.8) 
Eti 1 (c(fc - *) 2 « 2 + d(* - *)«* 2 + g) a fc _i Oi e fca2 , for fc > 1 

( for fc = 1,2 

if > = (4.9) 
I Ej^Eii^afc-iS-iate* * for A; > 2 

Comparing the coefficients of e kaz for each k, one has for k = 1: 

(a 4 - 6 a 2 +o) oi = 0. (4.10) 
Assuming ai 7^ (otherwise = for all k > 1 by induction), results in the values: 



"1-4 = ±y 2 , (4.11) 

which, in region 1 of FigfTl are of the form ai = —A + iu, 02 = «i, 03 = A + iu), 04 = S3. 
As our series solution (|4.5|) needs to converge for z > 0, we pick the two roots with 
negative real parts: 

00 

(z) = ^a k e kaiz + c.c. (4.12) 



it 



For A; > 1 one has: 



where the polynomial: 



p(ka l )a k = F[ k ' ) + F!?\ (4.13) 

p(k ai ) = (k ai )^ - b(k ai ) 2 + a, (4.14) 
is different from zero for k > 1, as may be seen from (|4,10p . Moreover, the terms into and 

(k) 

F 2 involves series coefficients aj, with j > A;, as can be easily observed from the expressions 



Therefore, the coefficients a&, k > 1 can be directly obtained from the equation (|4.13p as 
follows: 

ak = t Pka\, (4.15) 

where ip k , k > 1 are known functions given in Appendix IA.1I and depending on ol\ and the 
coefficients of the equation (|4.2|) . 
The first part of the homoclinic orbit corresponding to z > has thus been determined in 
terms of a±: 

00 

u + {z) = ai e aiz + fka k e kaiz + c.c. (4.16) 

k=2 
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We shall now construct the second part corresponding to z < 0. We remark that since the 
equation (|4.2p is reversible with respect to the transformation (|4.3|) . then for all z < we 
have that: 



oo 



u~(z) = b ie aiZ + Y, ip k b\e ka4Z + c.c. 



(4.17) 



k=2 

As we seek continuous solutions to (j4.2j) at z = 0, i.e. : 



oo oo 



U+(0) 



(4.18) 



it is sufficient to choose 61 



fc=2 fc=2 

ai and the relation 



00 



0-1 + Y = 



(4.19) 



k=2 

must hold. As the (pi~ are uniquely determined in terms of a\ and the coefficients of the 
equation (14. 2p . the relation (I4.19D allows one to compute a\ (in practice this must be done 
numerically at fixed values of the coefficients of the equation (|4.2p ). 

We have thus determined that the equation (|4.2[) has a homoclinic orbit to the origin of the 
form: 



{ u~(z) z < 

However, the converegence of the above series is still unclear. So as to not get mired in 
details and lose the main thread of our development, we defer proof of convergence to Section 
6, assuming its validity for the present. Note also that the above solution in (I4.20p represents 
homoclinic orbits at ALL (6, a) points in Region 1 of Figure 1. 

In Figj2]-(a) we fix a = 0.8 and b = 1.5 in region 1 and the other parameters are c = 
0.2, d = 0.1, g = 0.05 and h = 0.02. The series coefficient a\ has been chosen as the root 
ai = 40.4440 — 14.2061i of (I4.19p . In Figj2]-(b) all the equation parameters are chosen as in 
Figj2j The series coefficient a\ has been chosen as the root a\ = 41.6756 — 0.9564i of (|4.19p . 
In both figures the solution u is shown with respect to the variable x, for fixed t. 

In Figj3]-(a) we fix a = 0.8 and b = —1.5 in region 1 and the other parameters are as in 
Figj2j The series coefficient oq has been chosen as the root a\ = 14.814 — 41.6433i of (|4.19p . 
In Figl3]-(b) all the equation parameters are chosen as in Figl2j The series coefficient a\ has 
been chosen as the root a\ = —12.0892 — 44.4016i of (|4. 19[) . In both figures the solution u is 
again shown with respect to the variable x at fixed t. 

Notice that the jagged part of the pulses in Figj3] becomes smooth by changing the strength 
of the nonlinearities, as shown in Fig. HI where a and b are still the same as in Fig|3l but the 
other parameters are chosen as c = 0.7, d = 0.6, g = 0.8 and h = 0.4. 

In order to demonstrate the traveling- wave nature of our series solutions, Figure [5] plots the 
solutions at constant time t = 0, . . . ,4. In Figj5}-(a) we fix a = 0.8 and b = 1.5 in region 1 and 
the other parameters are c = 0.3, d = 0.1, g = 0.05, h = 0.02. The series coefficient a\ has 
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Figure 2: Symmetric home-clinic solutions. For the parameter choice see the text. 
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Figure 3: Symmetric homoclinic solutions: the pulses show a jagged part. For the parameter 
choice see the text. 



been chosen as the root a\ = -18.4550 + 15.8744i of (j4~T9l) . In FigEJ-(b) a = 3 and b = -1.5 
and the series coefficient a\ = 12.5579 + 36.8257i, the other parameters are as in Figj5]-(a). 
In addition to the convergence of our multi-infinite series solutions which, as mentioned 
earlier, is detailed in Section 6, we may also check our results versus earlier results obtained 
by other methods, for instance some derived numerically by homoclinic shooting [19]. Scaling 
the governing equation in [19], it is easy to see that only Figures 3 and 22 in that paper lie 
in Region 1 of our Figure 1. The homoclinic orbit given by our series solutions ()4.20p for 
parameters corresponding to those two figures in [19] are shown in Figj6] below, and they 
agree well with the numerical solutions in the earlier paper. 
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Figure 4: Symmetric homoclinic solutions: the jagged part of the pulses becomes smooth 
choosing stronger nonlinearities. 
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Figure 5: Traveling wave solution 



5 Non-reversible case: asymmetric homoclinic solutions 

Since the series solutions in the previous section agree closely with numerical results for the 
reversible case, we now include other nonlinear terms in (|2.ip which break the reversibility 
of the equation: 

/ // /// // try O n III II I O / 

f(u,u ,u ,u ) = cuu + du + gu + hu + puu + quu + ruu + su u . (5-1) 
The resulting equation is: 

mi ii n in 2 t 3 /;/ " ' S ; / \ 

u —bu + au = cuu + du + gu + hu + puu + quu + ruu + su u . (5-2) 

Let us build Shilnikov-type homoclinic orbits to the origin in region 1 of Figure 1, as in the 
previous section. Suppose that, for z > 0, the solution has the form (|4.5p . 
Let us substitute the series (|4.5p into the equation (|5.2p and it reduces to: 
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Figure 6: Comparison of our series solutions with the numerical solutions in [19]. In (a) the 
parameters are a = 7.5, b, d = — 3.75,c = —7.5,g = — 11.25, h = 0. In (b) the parameters are 
a = 15,6 = 3.75,c = —7.5, d = -3.75,5 = -11.25, /i = 0. 



with: 



where: 



J2 {(ka) 4 - b(kaf + a) a k e kaz = F 1 + F 2 + F 3 



k=i 



k=2 



F 2 = J^F 



(fc) 

2 > 



k=3 



f 3 = £k 



(k) 



fe=4 



(5.3) 



(5.4) 



F, 



(fc) 



F 



(*) 



£i=i - i ) 3 " 3 + (k- i) 2 (c + iga)a 2 + 



r(A; — i)a + d{k — i)ia 2 + g) a k - 



di e 



kaz 







Ej=2 ELl h a k-3 a 3-l a l e 

o 

Ei=s T!r=2 ECi 1 s a k-i a j-i a i a -i-i ekaz 



kaz 



for k = 1 
for A; > 1 

for k = 1,2 

for /c > 1 
for fc = 1,2,3 

for A; > 1 



(5.5) 



(5.6) 



(5.7) 



Comparing the coefficients of e kaz for each k, one has for k = 1 the equation (|4.10p . Once 
again we assume ai 7^ (otherwise = by induction) for all k > 1, resulting in the four 
complex roots fj4. 111) . 

As our series solution (|4.5p needs to converge for z > 0, we pick the two roots with 
negative real parts and the form of the solution for z > is still (|4.12p . Moreover, for k > 1 
one obtains the form of the coefficients a& as in (|4.15p . where 9?^, k > 1 are known functions 
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depending on ot\ and the coefficients of the equation (|5.2[) . The details are different now and 
are given in Appendix IA.2I 
The first part of the homoclinic orbit corresponding to z > has thus been determined in 
terms of ax, formally as in fj4. 161) . 

We shall now construct the second part corresponding to z < 0. Since the equation is 
not reversible we do not have any symmetry property for the solution, therefore impose the 
following solution form: 



oo 

u~(z) = J2he kaz + c.c. (5.8) 
fe=l 



where 9te(a) > 0, because the solution u (z) needs to converge for z < 0. Working as for 
z > 0, we obtain for k = 1 the following equation: 

(a 4 -ba 2 + a)h = 0. (5.9) 

Assuming bx ^ (otherwise bk = for all k > 1 by induction) , we find the four roots in (14. lip 
and, for z > 0, we choose a = txj. Moreover, for k > 1 one obtains the following equation: 



p{ka 4 )b k = F[ k) + F 2 (fc) + F 3 (fc) , (5.10) 



where the polynomial 



p(ka 4 ) = (kaf - b(kaf + a, (5.11) 
is non-zero for k > 1 and the quantities , and are formally the same as in 



(|5.5p - (15.7p once the coefficients at are substituted by bk- 
Therefore, the series coefficients can be easily obtained as follows: 

b k = i> k b\, (5.12) 

where Vfc are given in terms of &i and the coefficients of equations (|5.2p . Notice that in 
this case there is no relation between <p k and ^ k as in the reversible case (see the details in 
Appendix E2J). 
As we want to construct a solution of the following form: 

(u + (z) z>0 

z = (5.13) 

u~(z) z < 

which is continuous at z = 0, we impose: 



oo 



u+(0) = ai + J^ fc a^ = (5.14) 

k=2 

oo 

u-(0) = bx + J2^kb k x=0. (5.15) 

k=2 

Hence, we choose ax and bx as the nontrivial solutions of the following polynomial equations: 
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ax + ^<p k a k x = 0, (5.16) 

k=2 

oo 

&i + X>*&? = 0. (5.17) 

fe=2 

In practice the equations (|5.16|) - (|5.17p are numerically solved at fixed values of the coefficients 
of the equation (|5.2p . These solutions are not unique. In the following we show some numerical 
simulations of the homoclinic orbit (|5.13|) for different values of the parameters of the equation 

Q. 

In FigJT] the parameters a = 0.6 and b = 1.5 are chosen in Region 1 of Fig[TJ The other 
parameters are fixed as p = 0.1, q = 0.3, c = 0.2, r = 0.1, d = 0.3, g = 0.5, h = 0.2 and 
s = 0.1. In FigJT]- (a) the solutions of equations (|5.16p and (|5.17p are respectively chosen 
as ai = 2.3840 - 8.9933i and b x = 11.8609 + 3.8370f; in FigEj-(b) they are chosen as ax = 
-3.1533 - 8.6762i and b\ = -10.4363 - 6.2990i. 




Figure 7: Asymmetric homoclinic solutions. For the parameter choice see the text. 

In FigJH] the parameters a = 0.5 and 6 = 1 are chosen into the region 1 of FigJT] The 
other parameters are fixed as p = 0.01, q = 0.3, c = 0.02, r = 0.1, d = 0.3, g = 0.05, 
h = 0.2 and s = 0.1. In FigJ7]-(a) the solutions of equations (|5.16p and (|5.17p are respectively 
chosen as a\ = —7.6505 — 6.8582i and b\ = 4.7888 + 6.1539i; in FigfSl-(b) they are chosen as 
a x = 10.7954 + 1.8878* and h = 7.4002 - 2.35401 

In FigGHthe parameters a = 0.9 and b = —1.2 are chosen in Region 1 of FigJH The other 
parameters are fixed as p = 0.01, q = 0.02, c = 0.7, r = 0.1, d = 0.6, g = 0.8, h = 0.5 and 
s = 0.03. In Fig{J]-(a) the solutions of equations (|5.16p and (|5.17p are respectively chosen 
as ai = -8.5158 - 4.7092i and b x = -7.4718 - 7.3032i; in FigEJ(b) they are chosen as 
a x = -8.5158 - 4.7092i and b x = -6.7027 - 11.8976*. 

Further examples are shown in FigllOl In particular, in Fi^l01(a) the parameters are chosen 
as a = 0.9 and b = —1.0 in region 1 of FigHJ and p = 0.7, q = 0.6, c = 0.8, r = 0.5, d = 0.01, 
g = -0.02, h = 0.1 and s = -0.03. The roots of equations (|5TT6) - (|5TTjl are ax = -6.9709 
and b\ = —7.2113. In FietlOKb) the parameters are chosen as a = 0.9 and b = —1.6 in region 
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Figure 9: Asymmetric homoclinic solutions. For the parameter choice see the text. 

1 represented in FigHJ and p = 0.8, q = 0.6, c = 0.85, r = 0.7, d = 0.01, g = -0.02, h = 0.1 
and s = -0.03. The roots of equations (|5TT6>([5TTT|) are oi = -6.9709 and b x = -7.2113. 
In order to demonstrate the traveling-wave nature of our series solutions as in the previous 
section, Figure [TT1 plots the solutions at constant time t = 0, . . . ,4. Notice that we are only 
plotting the right solutions u + to avoid mess into the figure. The parameters are chosen 
as c = 0.3, d = 0.1,5 = 0.05, h = 0.02, p = 0.1, q = 0.2, r = 0.1, s = 0; in FigEKa) we 
fix a = 0.8 and b = 1.5 in region 1 and the series coefficient a\ = —16.0853 — 18.001H, in 
FigHEKb) a = 0.8 and b = -1.5 and a x = 4.8877 + 29.97081 

6 Proof of convergence 

We now return to the part of our derivation which had been deferred, viz. the proof that 
the series solutions obtained for the homoclinic orbits in the two previous sections are indeed 
convergent. In particular, as part of this proof, we will demonstrate the feature which was 
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Figure 10: Asymmetric homoclinic solutions. For the parameter choice see the text. 
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Figure 11: Traveling wave solution 

observed numerically earlier, i.e. that the later coefficients in our multi-infinite-series solutions 
drop off rapidly. 

As discussed in the previous two sections (see (|4.15p and f|5. 12j) respectively), one may 
inductively obtain the k-th coefficient in the series in the form 

a k = ip k a\. (6.1) 
More specifically (here we report the expressions into the symmetric case), for k = 2: 

{{c + d)a\ + g) a\ 



For k = 3: 



«3 



P(2 ai ) ' 

((5c + Ad)a\ + 2g) a\a 2 + h a\ 
p(3ai) 



(6.2) 



(6.3) 



15 



For k = 4: 



For k = 5: 



(10c + 6d)al + 2g) ai a 3 + (4(c + d)af + 5) a\ + a 2 , e 

a4 = Ri^) (6 ' 4) 



((17c + 8d)a? + 2g) a\a^ + ((13c + 12cf)a? + 2g) 0203 + 2>ha\a^ + 3/iaia2 

as = R^) (6 - 5) 

And so on for A; > 5: 

»* - <6 - 6) 

(A:) (k) i- 

where and -F 2 can be inductively obtained as the product of af with known functions 
of a\ and the parameters of the equation. 
Following some standard steps in |20j-|21j, it is straightforward to obtain the bound 



<rM,0- (fc+1 Vi| fc , (6-7) 



where r(c,d,£) is a constant dependent on the parameters c and d of our ODE, as well as 
Euler's constant £. Note that this validates the steep fall-off in the magnitude of the late 
terms of our series solutions which was numerically observed in Sections 4 and 5. 
For typical system parameters considered in Sections 4 and 5, the above bound may be 
expressed (evaluated) as: 



\a,h 



< W-^\ ai \ k ,k > 4. (6.8) 



Hence, |afc|,fc > 4 are bounded by some constant L > 0. Using this, we may bound the 
series of absolute values associated to the two halves (|4.5j) and (14. 17ft of our homoclinic orbit 
solutions by L times a convergent exponential (geometric) series. Hence our series solutions 
are absolutely convergent, and hence convergent, by the Comparison Test. 

As a final comment, the fixed point (b, a) = (0,0) of our ODE system is, as detailed earlier 
in Section 3, a saddle- focus in Region 1 of Figure 1. Hence we may conclude that our series 
solutions (14. 5b and (|4.17p represent the two halves of the homoclinic orbits for all (6, o) values 
in that region. 



7 Conclusions 

In summary, we have derived convergent, multi-infinite, series solutions for the homoclinic 
orbits of our canonical fourth-order ODE system in Region 1 of Figure 1, for both reversible 
and non-reversible cases. If the ODE is a traveling-wave reduction of some PDE, these 
solutions would correspond to genuine or localized pulses of the PDE. Given the many systems 
to which our canonical ODE pertains, our results are clearly of topical interest. In addition, 
the results for the non-reversible case go beyond those that have been typically considered in 
most analyses conducted within bifurcation-theoretic settings. 
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As discussed in Section 3, the homoclinic solutions in Region 3 of Figure 1 are known to be 
homoclinic or asymptotic to exponentially small periodic orbits (corresponding to delocalized 
solitary waves with exponentially small tails of any corresponding PDE(s)). Considering such 
solutions is clearly a natural next step for future work along the lines followed here for Region 
1 of Figure 1. One possible approach might be within the framework of |2J, treating our series 
solutions here as the unperturbed solutions uq in that work, and perturbing around them to 
derive the tail amplitudes of delocalized solitary waves in Region 3 of our Figure 1. Setting 
the tail amplitudes to zero would then yield the isolated curves in (6, a) parameter space 
(so far only known numerically) on which these delocalized waves become genuine embedded 
solitons [12]. However, the steps in the calculations would now involve a virtuoso use of 
multi-infinite series, and we defer this to a future paper. 

In another forthcoming paper the solutions and results here will be juxtaposed with, and 
placed in the framework of, a comprehensive categorization of ALL possible traveling-wave 
solutions, both smooth and non-smooth, for our governing canonical ODE. 

A Details of computations 

A.l Details for the symmetric case 
A. 2 Details for the general case 

For k = 2 the equation (|5.10j) reduces to: 

p(2ai)a 2 = ((p + q)a\ + (c + d)a\ + rai + g) of, (A.l) 
which gives a 2 in terms of a\, in the form given in (|4.15j) . where: 

(p + q)a\ + (c + d)a\ + roti + g 

<P2 = 77. — \ • (A.2) 

p{2ai) 

Analogously, for k = 3 the equation (15.10j) reduces to: 

p(3ai)a 3 = ((9p + 6g)af + (5c + M)a\ + 3rai + 2g) a±a 2 + ha\ . (A.3) 

Taking into account the expression for a 2 in terms of a\ obtained via ()A,6p . the equation 
(|A.3P gives 03 in terms of ai, in the form given in (|4.15p . where: 

9?3 = -77^ \3(3p 2 + 5pq + 2q 2 )at + (c(Up + Uq) + d(Up + Wq))al + 
p{3ai) 

(3r(4p + 3q) + c(5c + 9d) + id 2 + Wh)af + (g(llp + 8q) + (A.4) 
r(8c + 7d))a\ + {g(7c + 6d) - 4bh + 3r 2 )a 2 + hgron + 2g 2 + ah] 



And so on for k > 3: 



« - ' k^T ■ (A ' 5) 

(k) (k) (k) i 

where Ff', Ff J and Ff J can be inductively obtained as the product of with a known 
functions of a± and the parameters of the equation (|4.2p . 
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To compute the solution for z < the procedure is perfectly the same and the expressions 
of the coefficients are formally the same of a/% once substituted a.\ with a 4 = —a\. For 
example, at k = 2 the equation one obtains: 

p(2a 4 )b 2 = ((p + q)a\ + (c + d)a\ + ra 4 + g) b\, (A. 6) 

which gives b 2 in terms of b\, in the form given in (|5.12p . where: 

(p + q)a\ + (c + d)a\ + ra 4 + g , . . 

^ 2 = • (A - 7) 

These expressions involve both odd and even powers of a 4 , therefore there is no longer 
any symmetry linking the functions tpk and ^ (in the reversible case they were perfectly 
coincident (ff. = tpk, k > 1, because they involved only even powers of a± and 04). 
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